Constraints on power spectrum of density fluctuations from PBH evaporations 
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We calculate neutrino and photon energy spectra in extragalactic space from evaporation of 
primordial black holes, assuming that the power spectrum of primordial density fluctuations has a 
strong bump in the region of small scales. The constraints on the parameters of this bump based 
on neutrino and photon cosmic background data are obtained. 



I. INTRODUCTION 



o 
o . 

It is well known that for a sufficient production of primordial black holes (PBHs) in early Universe the spectrum 
of density perturbations set down by inflation must be "blue", i.e., it must have more power on small scales. This 
implies that the spectral index of the scalar perturbations must be larger than 1, in strong contradiction with the 
latest WMAP results [1]. In particular, inflationary models of hybrid type, in which the inflaton is trapped in a local 
minimum of the potential and which predict blue spectra seem to be excluded as a possible source of PBHs. 

As an alternative one can consider the models in which the power on small scales is boosted by means of a bump in 
the power spectrum of primordial fluctuations, as suggested, e.g., in [2]. Such a bump is predicted in many scenarios 
I , of two-step inflation (see, e.g., [3,4]) and can, in principle, exist even in single-field inflationary models (if the potential 
J> ■ has a special form leading to the period of fast-roll at the end of inflation [5] . 

Q\ \ In the present work we reconsider the problem of constraining the power spectrum of primordial fluctuations 
calculating the process of the formation of PBHs having small masses (10 11 — 10 15 g). Products of evaporation of 
these PBHs contribute, in particular, to extragalactic photon and neutrino diffuse backgrounds (which are measured 



experimentally) . 



II. PBH MASS SPECTRUM 



O 

Using the Press-Schechter formalism [6], the differential mass spectrum of primordial black holes can be written in 
O; the form [7,8]: 
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One must note, before the explanation of the notations used in Eq. (2.1), that in the problem studied here it is 
convenient to specify the spectrum of primordial fluctuations at a fixed time rather than at horizon crossing. It is so 
because we assume that, at the end of inflation, reheating and subsequent formation of the perturbation spectrum 
occur almost instantaneously. The connection between density contrasts at a fixed time and at horizon crossing 
(neglecting, at the moment, a dependence of the gravitational potential on the time) is quite simple [9], 
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Correspondingly, the mean square deviation (after smoothing the density field on a given mass M) at horizon crossing, 
<7h(M), is given by the formula 

M\4/3 f / k \4, ,..dk /Af\4/3 



Here, Afj is the horizon mass at the end of inflation, ai and Hi are the scale factor and Hubble parameter at the end 
of inflation, R is the comoving smoothing scale, R = 1/kn, connected with M by the expression 
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FIG. 1. Horizon crossing amplitude 8u{k) (dashed line) and smoothed dispersion an(kfi) (solid line). Spectra are shown for 
S = 3,6% = 0.06, k = 2.75 ■ 10 16 Mpc' 1 . 



The notations used for the other functions in the integral (2.3) are standard: 5u{k) is the horizon crossing amplitude 
for primordial perturbations of the density contrast, W(kR) and T(k) are window and transfer functions. 

The variable 5^ in Eq. (2.1) is a value of the (smoothed) density contrast at horizon crossing, pi is the energy 
density at the end of inflation. At last, the function f(M, <5|f ) connects the values of the smoothing mass M, density 
contrast <5|f and PBH mass Mbh, 



M BH = f(M,5%;M l ). 



(2.5) 



The concrete form of this function depends on the features of the gravitational collapse leading to PBH formation. 
It is assumed that the dependence of Sff(k) on k at small scales can be parameterized in the form 



\gS H (k) = A+(lgd° H - A)exp 
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Here, the value of A is known from observations at large scales, 5%, ko and E are free parameters, the values of which 
should be constrained. 

The calculation of <jh{M) was performed using the gaussian window function, 
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and the transfer function, derived in the cosmological perturbation theory, 
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r is the conformal time. In the following we will ignore the time dependence of the transfer function, putting 

1 
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(2-10) 



(because the most abundant PBH production takes place at the time when the scale k^ 1 corresponding to the 
maximum of Sn(k) crosses horizon). 

In the model of spherically-symmetric collapse one obtains the expression 



M BH = (6%) 1 ' 2 M h , 



(2.11) 
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FIG. 2. PBH mass spectra calculated for the models of critical collapse and Carr-Hawking collapse. Here, we used the 
following set of parameters: £ = 3, 5% = 0.06, ko = 2.75 • 10 Mpc~ , Trh = 10 GeV. The parameters for the critical collapse 
were taken as follows: S c = 0.45, j c = 0.36, k c = 4. 



where M h is the fluctuation (smoothing) mass at the moment of horizon crossing. According to Carr and Hawking 
[10], 1/3 <6% < 1. 
The connection between M and the horizon mass M/j is [7] 

M h = M- /3 M 2/3 . (2.12) 

From Eqs. (2.11), (2.12) one has the expression for the function f(M,S^; Mi) for the Carr-Hawking collapse: 

M BH = f(M, 5% ; Mi) = (^) 1 / 2 M 2 / 3 M^ /3 . (2.13) 

In the picture of the critical collapse [11] the corresponding function is 

f(M, 5% ; Mi) = k c {5" - 6 c )^M 2 / 3 Ml /3 , (2.14) 

where 5 C , j c and k c are model parameters. In this work we will accept the following set of parameters, which is in 
agreement with recent calculations [12]: 

5 C = 0.45, 7c = 0.36, k c = 4. (2.15) 

Some results of PBH mass spectrum calculations are shown in Fig. 1 and 2. In Fig. 1, the smoothed dispersion is 
drawn as a function of the inverse comoving smoothing scale, kfi — l/R. The connection between kfi and horizon 
mass Mh is determined from Eqs. (2.4) and (2.12) and is 



Fig. 2 shows that the resulting PBH mass spectrum strongly depends on the model of the gravitational collapse 
(at the same values of the fluctuation spectrum parameters and reheating temperature T RH ). 
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FIG. 3. Red shift distribution of the integrand F(E, z). Left panel is for neutrinos, right panel - for photons (absorption of 
7-rays at z > 700 is not shown in this figure). PBH mass spectra shown in Fig. 2 were used, E = lGeV for both panels. 



III. NEUTRINO AND PHOTON SPECTRA FROM PBHS EVAPORATIONS 

Evolution of a PBH mass spectrum due to the evaporation leads to the approximate expression for this spectrum 
at any moment of time: 

2 

where a accounts for the degrees of freedom of evaporated particles and, strictly speaking, is a function of a running 
value of the PBH mass m. In all our numerical calculations we use the approximation 

a = const = a(M%% x ), (3.2) 

where Mg^f is the value of Mbh in the initial mass spectrum corresponding to a maximum of this spectrum. Special 
study shows that errors connected with such an approximation are rather small. 

The expression for an extragalactic differential energy spectrum of neutrinos or photons (the total contribution of 
all black holes) integrated over time is [7] 



= J F(E, z)dlg(z + l). 



(3-3) 



In this formula, dj, a, and (Xq (XY& cosmic sa lie factors at ti, t and at present time, respectively, and ip(E,m) is 
a total instantaneous spectrum of the radiation (neutrinos or photons) from the evaporation of an individual black 
hole. The exponential factor in Eq. (3.3) takes into account an absorption of the radiation during its propagation in 
space. The processes of neutrino absorption are considered, in a given context, in [7]. In the last line of Eq. (3.3) we 
changed the variable t on z using the flat cosmological model with Oa / for which 
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dz H (l + z) 
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with Q r = (2.4 • 10 4 /i 2 )-\ h = 0.7, Q M = 0.3, n A = l-Q M - Q r . 

The instantaneous spectra of neutrinos and photons from PBH evaporations were calculated using the photosphere 
model elaborated in works [13,14]. The decoupling temperature Tf for photons was taken equal to 120 MeV, and the 
temperature of the ncutrinosphcre, T v , is equal to 100 GeV, as in [14]. The neutrino emission from "cold" PBHs, 
those with Th < 200 GeV, was calculated using the standard Hawking formula. The contribution to the black hole 
neutrino spectra from decays of pions and muons was, at this stage, not taken into account. 

The photosphere model of black hole evaporation used here predicts very steep time-integrated spectra of photons 
and neutrino (E~ 4 ), and in the low energy region the absolute spectrum values (the number of particles) are very 
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FIG. 4. Constraints on the horizon-crossing amplitude <5h(A;) that were obtained in this work. M° is the horizon mass 
at the moment when fluctuation with comoving wave number ko enters horizon. Constraints following from both neutrino 
and gamma ray experiments are shown. Dashed lines correspond to the model of critical collapse with the parameters 
S c — 0.45, 7c = 0.36, k c — 4, solid lines represent the results obtained using the standard collapse picture. 

high (comparing with the models without photosphere) . The average energy of photons and neutrinos as a function 
of PBH lifetime is proportional to the square root of black hole temperature (E <~ \/Th rather than E ~ Tjj). 

On Fig. 3 the red shift distributions of the differential energy spectra S(E) are shown for the neutrino or photon 
energy 1 GeV. One can see that the distributions are more wide in the case of the critical collapse (because the PBH 
mass spectrum in this case has a long "tail" of small masses). The absorption factor e~ T is sufficient: for photons 
Zmax ~ 700 [15] and for neutrinos it is about 10 6 [7]. 

IV. RESULTS AND DISCUSSIONS 

The calculation of the constraints uses two basic observational facts. 

1. The differential energy spectrum of the extragalactic photon background at energy E 1 = lOMeV is [16] 

- lO^GeV^cm^s^sr- 1 . (4.1) 

2. According to the data of Super-Kamiokande experiment [17], the electron antineutrino background flux in space 
is constrained by the inequality 

^(E Pe > 19.3GW) < 1.2 cm- 2 .?' 1 . (4.2) 

The resulting constraints are shown on Fig. 4. The forbidden values of the most important parameter, fo(fco), he 
inside the bordered regions drawn on the figure. The value of the parameter £ characterizing the width of the gaussian 
distribution in Eq. (2.6) was fixed throughout all calculations (£ = 3). The constraints are given as a function of 
the horizon mass corresponding to the moment of time when the comoving wave number ko enters horizon. It follows 
from these results that the constraints are stronger in the case of the standard Carr-Hawking collapse. 

The second important result is that the constraints based on neutrino emission of PBHs are comparable with those 
following from photon emission. At the region of small horizon masses (large ko) where the large red shifts are 
sufficient, the constraints from neutrino emission are stronger. 

One should add, at the end, one comment. In this work we calculated the constraints using the rather high value 
of Trh (= 10 10 GcV). The corresponding horizon mass at the end of inflation, Mi, is equal to ~ 10 11 g. The interval 
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of studied here values of the kg parameter corresponds to PBH production in super-horizon perturbations. At low 
Trh the value of M, will be high and the parameter kg can be chosen such that black holes will be produced from 
sub- horizon perturbations (the possibility of PBH production from such perturbations is studied in [18]). In principle, 
constraints on the power spectrum due to PBH evaporations can be obtained in this case also, using the hypothesis 
of critical collapse. 



[1] D. N. Spergel et al, arXiv:astro-ph/0603449. 

[2] P. Ivanov, P. Naselsky and I. Novikov, Phys. Rev. D 50 (1994) 7173. 

L. Randall, M. Soljacic and A. H. Guth, Nucl. Phys. B 472 (1996) 377. 

J. Garcia-Bellido, A. D. Linde and D. Wands, Phys. Rev. D 54 (1996) 6040. 

T. Bringmann, C. Kiefer and D. Polarski, Phys. Rev. D 65 (2002) 024008. 
[3] M. Yamaguchi, Phys. Rev. D 64, 063503 (2001). 

[4] M. Kawasaki, T. Takayama, M. Yamaguchi and J. Yokoyama, Phys. Rev. D 74, 043525 (2006). 

[5] S. Chongchitnan and G. Efstathiou, arXiv:astro-ph/0611818. 

[6] W. H. Press and P. Schechter, Astrophys. J. 187 (1974) 425. 

[7] E. V. Bugaev and K. V. Konishchev, Phys. Rev. D 65 (2002) 123005. 

[8] E. V. Bugaev and K. V. Konishchev, Phys. Rev. D 66 (2002) 084004. 

[9] E.W. Kolb and M.S. Turner, The Early Universe (Addison- Wesley, New York, 1990). 
[10] B. J. Carr and S. W. Hawking, Mon. Not. R. Astron. Soc. 168, 399 (1974). 
[11] J. C. Niemeyer and K. Jedamzik, Phys. Rev. Lett. 80 (1998) 5481. 

J. C. Niemeyer and K. Jedamzik, Phys. Rev. D 59 (1999) 124013. 
[12] I. Musco, J. C. Miller and L. Rezzolla, Class. Quant. Grav. 22, 1405 (2005). 
[13] R. G. Daghigh and J. I. Kapusta, Phys. Rev. D 65, 064028 (2002). 
[14] R. G. Daghigh and J. I. Kapusta, Phys. Rev. D 67, 044006 (2003). 
[15] A. A. Zdziarski and R. Svensson, Astrophys. J. 344, 551 (1989). 
[16] A. W. Strong, I. V. Moskalenko and O. Reimer, arXiv:astro-ph/0506359. 

[17] M. Malek et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 90 (2003) 061101 [arXiv:hep-ex/0209028]. 
[18] D. H. Lyth, K. A. Malik, M. Sasaki and I. Zaballa, JCAP 0601 (2006) 011. 
I. Zaballa, A. M. Green, K. A. Malik and M. Sasaki, arXiv:astro-ph/0612379. 



G 



